Boundary Hamiltonians for topological models and Projected Entangled Pair States 
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In the ground state of a gapped spin system, the physical properties of a region of the bulk are 
encoded in a Hamiltonian that lives on its boundary. We study the structure of such boundary 
Hamiltonians for topological models in the framework of Projected Entangled Pair States (PEPS). 
We find that the boundary Hamiltonian decomposes into two parts: A universal part which is entirely 
non-local but independent of microscopic details, encoding the nature of the topological phase, and 
a non-universal part which is local and depends on microscopic details only. A topological phase 
transition is reflected by a diverging interaction length of the non-universal part. 



The study of strongly correlated quantum systems is 
of central interest in modern condensed matter physics 
due to the exciting properties exhibited by those sys- 
tems, in particular unconventional phases with topolog- 
ical order. In order to identify topological order in such 
systems, topological entropies [J, Q have been applied 
successfully. To obtain more information about the en- 
tanglement of the system than contained in the entropy, 
the entanglement spectrum (ES) - i.e., the spectrum of 
the reduced density operator of a region - has been stud- 
ied, and it has been realized that for certain systems, the 
low-energy part of the ES resembles the spectrum of the 
thermal state of a one-dimensional (ID) local boundary 
Hamiltonian which can be associated to the boundary of 
the region studied, and which seems to be related to the 
edge excitations present in the model 0-111 ■ While this 
relation between bulk entanglement spectrum, boundary 
Hamiltonian, and edge excitations can be made rigorous 
in some cases 0, 0| j in rnost cases the boundary Hamilto- 
nian is a posteriori inferred from the structure of the en- 
tanglement spectrum, and a general connection between 
entanglement spectrum and boundary still needs to be 
made. 

In Ref . , we made progress in that direction by prov- 
ing a rigorous connection between entanglement spec- 
trum and boundary using the framework of Projected 
Entangled Pair States (PEPS) which form the ap- 
propriate description of ground states of gapped lo- 
cal Hamiltonians both in conventional and topological 
phases [llj. We found that in trivial phases (without 
symmetry breaking or topological order), the boundary 
Hamiltonian is local, while it becomes non-local at phase 
transitions, allowing to identify phases by looking at the 
locality structure of the boundary Hamiltonian. We also 
found that for topological models such as the toric code 
or quantum double models 1^ which are renormaliza- 
tion group (RG) fixed points, the boundary Hamilto- 
nian becomes completely non-local, coupling all spins si- 
multaneously. However, the question remained open of 
how to identify topological phases away from RG fixed 



points from their boundary Hamiltonian, as well as how 
to detect phase transitions between topological and non- 
topological phases. 

In this paper, we establish a framework for studying 
the boundary Hamiltonians of systems with topologi- 
cal order in the framework of PEPS. We find that the 
boundary Hamiltonian of topological models splits into 
two parts. The universal part determines the phase of 
the system: It couples to a topological degree of free- 
dom at the boundary (such as the parity) and is thus 
completely non-local, it couples boundaries at arbitrary 
distance, and it depends only on the phase of the sys- 
tem and not on the exact point within the phase. On 
the other hand, the non-univeral part is local (generaliz- 
ing what happens for trivial phases), it does depend on 
the point within phase, and it vanishes under RG flows. 
At the phase transition, the non-universal part becomes 
non-local, thereby destroying the topological order con- 
tained in the universal part, again allowing to detect the 
phase transition at the boundary. 

Let us first introduce PEPS and explain how to use 
them to derive ID boundary theories from 2D models. 
For clarity, we restrict to square lattices on a cylin- 
der (with length Nh and circumference N^). A (trans- 
lational invariant) PEPS \Tp) = '^ii---ijv Nii • ■ • i *Jv) is 
described by a five-index tensor A^^^^g (Fig. [1^, with i 
called the physical and a, /3, 7, 6 the virtual indices), such 



that the coefficient c, 



is obtained by arranging ten- 




FIG. 1. Tensor networks (see text for details). 
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sors , • ■ • , on the cyUnder and contracting each 
virtual index with the corresponding index of the ad- 
jancent tensors, and putting some boundary conditions 
IX-l)i \xh) at the open virtual indices at the two ends 
(Fig. [It). As proven in Ref. for any PEPS the entan- 
glement spectrum of a half-cylinder (green cut in Fig.[TJ;) 
is equal to the spectrum of ct (x a/ct^ cr/j ^/a^ (with * the 
complex conjugate), where an is the state obtained at 
the virtual indices of the right half-cylinder by contract- 
ing the physical indices with the adjoint, cf. Fig. [TJi,f 
(with xr = \xr){xr\, and correspondingly for ctl). In 
turn, this allows to construct a boundary Hamiltonian 
H = — log a which is associated to the virtual degrees of 
freedom at the boundary, and which exactly reproduces 
the entanglement spectrum. As demonstrated in [§], the 
phase of a system can be detected by looking at the local- 
ity of its boundary Hamiltonian: It is short-ranged in the 
trivial phase; the interaction length diverges at the phase 
transition, giving rise to a diverging correlation length; 
and it is long-ranged in non-trivial phases. Naturally, 
boundary Hamiltonians with short-range interactions are 
of particular interest, as they can be related to physical 
models. In this paper, we will reconsider the boundary 
Hamiltonian for non-trivial phases, and show how a local 
boundary Hamiltonian can be recovered. 

Topological order in PEPS is accompanied by a gauge 
(i.e., virtual) symmetry of the tensor A, such as the in- 
variance under the representation of a (finite) symmetry 
group. Fig. [TId (more general symmetries are give n, 
e.g., by Hopf algebras [IJ] or tensor categories [l5|). 
For simplicity, we will focus on Z2 symmetry, i.e., A 
is invariant under Z^'^, for some unitary representation 
{^,Z} of Z2, but our findings easily generalize to any fi- 
nite group. The simplest instances of topological models 
are renormalization (RG) fixed points, such as the toric 
code [I4I for Z2. There, the tensor A* is isometric on 
the symmetric subspace [l3| . i.e., the transfer operator 
E [Fig. [TJi] equals the projector onto the symmetric sub- 
space, E = i(l®^ -I- Z®"'), when understood as a map 
from ket to bra indices. This form of E is stable under 
contraction with itself (reflecting the fact that we are at 
the RG fixed point), and thus, one column of the trans- 
fer operator (Fig. [Ih) equals T = i(l®27V„ ^ ^»2Ar„-)^ 
Again, we have that = T as we are at the RG fixed 
point, and thus, the two virtual boundaries of Fig. [T^ 
are jointly described by aLR oc l^^JVt. _|_ ^mN^ ^ which 
corresponds to a completely non-local boundary Hamilto- 
nian H^fj — PZ®'^^^ coupling the two boundaries, with 
/? = 00. Imposing a boundary condition on one side 
(Fig. [If) leads to (Tr/l oc al^^" -I- bZ®^\ where a and 
h depend on the corresponding boundary condition, and 
thus yields an boundary Hamiltonian H = /SZ"^^' where 
(3 depends on the boundary conditions, cf. Ref. [9[. 

To gain a better understanding of the structure of this 
boundary, it will be helpful to rewrite it as follows: Define 



the projectors onto the even and odd parity subsectors, 
respectively. Then, we can rewrite aLR. = f oven ^ -Fovcn + 
-Podd ^ Podd, i.e., the boundary in Fig. [T^ is such that 
the total parity of the left and right boundary is even — 
the boundary can thus be understood as imposing a par- 
ity superselection r ule0 Correspondingly, by imposing 
a boundary condition XR on the right side (Fig. [If), we 
obtain a boimdary state an = PcvcnPcvcn + pfdd^odd, 
with = 1 - pf^j = tr[PcvcnXfl]- By also fixing 

a left boundary condition xl (with 75^011 and p^^jj de- 
fined accordingly), we find that the boundary Hamilto- 
nian is H = - logPovcn-Povcn ~ logPodd-Podd (with Pevcn = 

PcvcnPfvon' ^ud Podd = PoddP^dd)' Thus, the boundary 
Hamiltonian at the RG fixed point acts by changing the 
weight of the parity subsectors, with the weight given 
by the boundary conditions. (Note that the choice of the 
boundary x in the PEPS can be achieved by choosing ap- 
propriate boundary terms in the corresponding physical 
Hamiltonian.) 

We would now like to go beyond the RG fixed point and 
study how topological order is reflected in the boundary, 
and whether there is a way to detect topological phase 
transitions and the disappeareance of topological order. 
In order to study these phenomena, we will consider the 
toric code with string tension which is obtained by in- 
troducing a specific magnetic field. This model has been 
first discussed in [ist . where it was shown that it ex- 
hibits a topological phase transition by mapping it to 
the 2D classical Ising model; see also [19[. The ground 
state of the model has an exact PEPS representation 
which can be constructed from the tensor of the toric 
code (cf. Ref. [31, Eq. 7.5) by acting with an invertible 
map A which supresses |1) states (strings) by a weight 
A < 1. Here, A = 1 corresponds to the toric code, and 
A = to the vacuum; the phase transition occurs at 
A = l/Vl + \/2 « 0.644 Ol. (Note that since A is 
invertible, we can kick it back to the toric code Hamilto- 
nian to obtain an exact frustration-free Hamiltonian for 
the deformed model [20j). 

The presence of a gauge symmetry. Fig. [If), in the ten- 
sor is the only known way to achieve topological order in 
PEPS. Let us therefore try to understand what the pres- 
ence of such a symmetry implies for the structure of the 
boundary. For now, we will focus on a cylinder of fixed 
circumference Ny , where we will take its length Nh to in- 
finity. By blocking the tensors in a column, we obtain a 
description of the state on the cylinder as a Matrix Prod- 
uct State (MPS) [2I, 22| (i.e., a one-dimensional tensor 

T^^. T inher- 



network). characterized by a tensor T 



(1 



Z«5^") and Podd = i(l^ 



This parity can be understood as arising from fermions with 
internal degrees of freedom: E^g., the toric code has Z = 
(^_j^), i.e., spinless fermions [l^ . and the RVB state has 
Z = diag(l, 1, — 1), i.e., spinful fermions with a hardcore con- 
straint such as in the t- J- model lid llTll . 
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its the symmetry from the PEPS tensors, = ZT'^Z, 
with Z = Z^^^', i.e., the matrices are block-diagonal, 
with the two blocks corresponding to the even and odd 
parity subsector, respectively. It follows that the 
transfer operator T = (g) T', Fig. [T^, has positive 

maximal left (right) eigenvectors in both sectors, cr^^^^ 
((T^en) ^iid cTqjjj (cr^d)' ^^^^ Corresponding eigenvalues 
7e and 70. Using the Cauchy-Schwarz inequality, it is 
straightforward to see that the largest eigenvalue |7eo| of 
the off-diagonal block is bounded as |7eo| < ^/je'^o- This 
is, for a sufficiently long cylinder the boundary states 
ctl and cr/j, and thus a cx \f^'^R\f^-, will be block- 
diagonal (up to an exponentially small correction) for any 
boundary condition, a = PovonO'cvon ® Poddg odd- ( Here , 

and trcTovcn = trcodd = 1-) Again, povcn and podd de- 
pend on the boundary conditions x- Conversely, for a 
given Nh we can always choose boundary conditions x 
such that any such a is realized. 

As we have seen, the virtual boundary of a system with 
topological gauge symmetries has a block-diagonal struc- 
ture a = pevenO-even © PoddCTodd, depending on the parity 
of the boundary conditions x- How should we constuct a 
boundary Hamiltonian from this? A first idea is to choose 
a boundary condition with fixed parity (w.l.o.g. even), so 

that a cx CTeven, and write (Teven OC exp[— i?]. Since (Teven 

is only supported on the even parity sector, H must then 
be of the form H ~ /3topo^odd + H' , with /3topo = 00 and 
H' = — log CTcvcn (with the log taken on the even par- 
ity subspace). Unfortunately, numerical study quickly 
reveals that the H' defined this way is highly non-local 
both in the topological and the non-topological phase 
(see inset of Fig. [2]); in particular, it has strong iV„ — 1 
and iV„-body terms. It is in fact easy to see where this 
comes from: Since H' = H'Z'^^^ , any local term h in H' 
comes with a "twin" hZ'^^^ which is necessarily highly 
non-local! There is an ad hoc solution to that by not- 
ing that H' is not uniquely defined, but only its even 
parity projection PcvcnH' Peven = - log (Teven- ThuS, WC 
can declare the most local H' among all possible choices 
the "right" one (this corresponds to counting each pair 
of interactions h and hZ^^^' with the locality of h). This 
solution does indeed show the desired behavior: It gives 
a local Hamiltonian H' in the topological regime which 
becomes non-local at the phase transition. 

Yet, this ad hoc definition of H' might feel unsatisfac- 
tory. Fortunately, there is a more elegant solution. To 
this end, recall that for arbitrary boundary conditions, 
the boundary state is given by cr = PovonCovon+PoddCodd- 
If wc now consider H = — log cr, it follows that it is of 
the form 

H = i^topo + H' , (1) 

with 

-fftopo = - log Povon ^'ovon - logPodd ^'odd (2) 



and 

= - log f ^ . (3) 

\ CTodd/ 

We thus find that the Hamiltonian splits into two parts: 
The first part, i/topo, depends on the choice of bound- 
aries Xj and it controls the relative weight of the even 
and odd parity subspace and is therefore nonlocal — in 
fact, it is exactly the same as the one we found for the 
toric code earlier. The second part, H' , on the other 
hand, is independent of the boundary conditions %; it 
only depends on the fixed points of the transfer opera- 
tors in each parity sector and is thus uniquely defined. 
Since i?topo = for Povcn = Podd, H' = H equals the 
actual entanglement Hamiltonian for appropriate bound- 
ary conditions; moreover, any boundary condition with 
Pcvcn 7^ Podd will generally result in a less local H since 
Htopo OC Z®^" -I- const. The numerical results for the lo- 
cality of H' for the toric code with string tension are 
shown in Fig. [2j We find that H' is short-ranged in 
the topologically ordered regime, becomes long-ranged 
at the phase transition, and remains so in the entire non- 
topological regime. We have then used this finding to 
study the interpolation between the square lattice Toric 
Code and the RVB state on the kagome lattice described 
in [l^l. The result is shown in Fig. [3l The boundary 
Hamiltonian H' is local throughout, and no sign of a 
phase transition can be found, thus strongly supporting 
the existing evidence that the RVB state is a topological 




interaction range 



FIG. 2. Interaction range of the boundary Hamiltonian Q 
for the toric code with magnetic field, for = 12. A = 1 is 
the toric code, A = the trivial state, and A « 0.644 the phase 
transition (red). Blue (green) denote the topological (trivial) 
phase. The interaction range of the boundary Hamiltonian 
(as defined in Ref. [^) decays exponentially in the topological 
phase, while it is long-ranged in the trivial phase. (To ex- 
clude finite-size effects, only points up to = 6 should be 
considered for the decay.) The inset shows the Hamiltonian 
defined only from (Tevon, which has strong 12-body terms in 
both phases. The calculations (as well as the ones in Figs. [3] 
and [4]) have been carried out using exact column- wise con- 
traction (cf. d, [H, [13) and are thus exact. 
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FIG. 3. Interaction range for the Toric Code-kagome RVB in- 
terpolation [l^ for — 8 (large symbols) and A''^, = 6 (small 
symbols), (a) shows the strength of fc-body terms, while (b) 
and (c) analyze the decay of two specific two-body terms (see 
Ref. [131) with distance. It is found that all terms decay expo- 
nentially throughout the interpolation, showing the absence 
of a topological phase transition between the toric code and 
the RVB system. 



Z2 spin liquid. 

The decomposition (jlj of the entanglement Hamilto- 
nian into -fftopo and H' can be understood as follows: 
Hto-po consitutes the universal part of the boundary 
Hamiltonian — it controls the relative weight of the par- 
ity sectors and is therefore completely non-local, does not 
depend on the microscopic details of the phase (here. A), 
is invariant under renormalization, and captures long- 
range features such as the fact that the two boundaries 
are connected. H' , on the other hand, is non-universal — 
it depends on the microscopic details of the phase, and it 
vanishes under renormalization in the topological phase; 
on the other hand, it becomes long-ranged at the phase 
transition and thus destroys the topological order other- 
wise contained in i?topo- 

What happens if we consider a PEPS cylinder with 
two virtual boundaries. Fig. [T^, which is what we e.g. 
encounter if we want to study the ES on a torus? If 

'^L = Cr*R, we have that T°° = (Teven ® (^even + CTodd CTodd, 

which impUes that the boundary Hamiltonian is of the 

form 

H = Ptopo{Pcvcn(^Pcvcn + Podd®Podd) + H'(g)l + l(g)H' , 

where /3topo — 00, and H' again as in Eq. ([3]): The uni- 
versal topological part couples the two boundaries and 
imposes that they joinly have even parity, while the non- 
universal part is the same local Hamiltonian as before, 
acting independently on both boundaries. ( If cr^ 7^ ajt, 
H' is different on the left and on the right.) 

The form ([T]) of the boundary Hamiltonian has also 
consequences for the topological entropy S{pl) = S{a): 
We have that 

S{Pl) = HiPcvcn) +PcvcnS'((Tcvcn) + Podd<S'(crodd ) , (4) 

where H{p) ~ — plogp— (1 — p) log(l — p), i.e., depending 
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FIG. 4. Scaling analysis (cf. text), (a) Diagonal/off diagonal 
gap for the toric code with string tension, (b) Same for the 
interpolation between square lattice toric code and the RVB 
state on the kagome lattice [l^. (c) Splitting between even 
and odd block for the Toric Code-RVB interpolation, (d) Dif- 
ferent of the phase-corrected nearest-neighbor terms for the 
two boundary Hamiltonians -ff" and . In all cases, the in- 
terpolation parameter < A < 1 is increased in steps of 0.05. 
7e, 7o, and 7eo are the largest eigenvalue of the even, odd, 
and off-diagonal block of the transfer operator, respectively. 



on the choice of the boundary condition, the topological 
correction can be maximal (peven = 1 or podd = 1), zero 
(Povon = Podd = 5 ) oi' anywhere inbetween. 

Up to now, we have considered systems in the limit 
where Nh — >■ oo for finite N^; we exploited this in par- 
ticular when deriving the block-diagonal structure of the 
boundary. Does this limit correctly reproduce the physics 
of systems where we take the simultaneous limit in 
and Njjl In order to understand this, we need to study 
the properties of the transfer operator as a function of 

Firstly, in order to obtain a block-diagonal boundary 
operator, it is important that the off-diagonal blocks are 
supressed relative to the diagonal ones for Ny — > oo. 
Fig. 0^ shows this data for the toric code with string 
tension: In the topological phase (blue), the supression 
converges to a constant as required; in the trivial phase 
(green), on the other hand, the supression vanishes ex- 
ponentially in Njj, and our limit is no longer justified. 
Fig. shows the same data for the toric code-RVB in- 
terpolation, proving that the off-diagonal blocks are in- 
deed supressed throughout. 

Secondly, if the eigenvalues corresponding to the even 
and odd parity subsector are not equal, in principle 
only one of them contributes to the boundary state for 
Nh ^ oo. For the Toric Code with string tension, there is 
no such splitting due to symmetries; for the Toric Code 
- RVB interpolation. Fig. [4]; proves that this splitting 
vanishes exponentially in N^. (Note though that even 
with a finite splitting, all topological sectors will appear 
as ground states in the thermodynamic limit as they are 
not coupled by the Hamiltonian [23j). 
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Up to now, we have only considered one of the two 
ground states on the cyhnder which exist for a fixed 
boundary condition. The second ground state can be ob- 
tained by moving one excitation from one boundary to 
the other; in the PEPS representation, this is achieved 
by placing a string of Z's connecting the boundaries on 
the bonds This already suggests that the difference 
between the boundary Hamiltonians lies in the presence 
of a TT flux through the torus, which gives rise to a Berry 
phase: We expect that the boundary Hamiltonian H'^ 
with Z string is the same as the one without, H'^ , except 
for a phase e™ for any term which describes hopping 
across the location of the Z string. Fig. |3Ji shows the 
difference between the nearest neighbor terms in and 
after undoing the e" phase in H'" , which indeed van- 
ishes exponentially with N.^ in the topological phase; the 
same behavior is observed for longer-range terms. Taking 
superpositions of these two ground states can be under- 
stood as entangling the flux |0) with i?*; this increases 
(jll) (and thus decreases the topological correction) by the 
mixing entropy. 

To conclude, we have studied the boundary Hamilto- 
nian of topologically ordered models and found that it 
consists of two parts: A non-local universal part which 
encodes the topological features of the system, and a lo- 
cal non-universal part which depends on the details of 
the system within the topological phase. The topological 
order is destroyed when the non-universal part looses its 
locality. The structure of the universal part is given by 
the symmetry underlying the model, which can be inter- 
preted in terms of particles with conservation laws which 
are given by the topological model under consideration 
(such as a group for double models, of a general fusion 
rule for string nets). It is an extremely suggestive idea 
that such a decomposition should exist for any phase, i.e., 
that for any phase there exists a universal contribution of 
the boundary Hamiltonian such that the remaining non- 
universal part becomes local, and this very universal part 
would then label the phase of the sytem. 
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